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ABSTRACT 

Several  techniques  for  finding  approximate  solutions  to 
certain  classes  of  nonlinear  differential  equations  are  in- 
vestigated.  The  nonlinear  systems  evaluated  are  second  order 
with  quadratic  and  cubic  terms,  having  driving  functions  and 
initial  conditions.   Pipes'  technique  is  used  to  reduce  the 
nonlinear  differential  equation  to  a  system  of  linear  dif- 
ferential equations.   The  Brady-Baycura  technique  makes  use 
of  nonlinear  Laplace  transforms  to  obtain  the  solution.   The 
solutions  are  compared  to  the  well  known  Runge-Kutta  numerical 
method  solution  using  the  digital  computer.   Greater  accuracy 
was  found  using  the  Brady-Baycura  method,  but  the  simplicity 
of  Pipes'  method  makes  it  more  attractive  to  the  engineer. 
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I.   INTRODUCTION 

Although  numerical  methods  such  as  Runge-Kutta  have  been 
widely  used  to  solve  nonlinear  differential  equations  by  means 
of  digital  computers,  the  engineer  does  not  always  have  access 
to  a  computer.   Practical  techniques  to  assist  the  engineer  in 
obtaining  approximate  solutions  are  needed. 

Pipes  (Ref.  1)  developed  a  method  of  solving  nonlinear 
differential  equations  by  a  reversion  technique  quite  akin 
to  algebraic  series  reversion.   The  technique  uses  Laplace 
transforms  extensively  in  obtaining  the  approximate  solution. 
Baycura  and  Brady  (Refs.  2  and  3)  have  proposed  and  examined 
nonlinear  Laplace  transforms  to  obtain  solutions  to  nonlinear 
differential  equations.   Tou ,  Doetsch,  Pipes  and  Nowacki 
(Refs.  4,  5,  6  and  7)  have  proposed  an  iterative  method  using 
Laplace  transforms  in  each  iteration.   Flake  (Ref.  8)  has  pro- 
posed a  cumbersome  method  using  Volterra  series  and  special 
transforms. 

The  techniques  were  investigated  and  applied  to  driven 
second  order  systems  with  initial  conditions.   The  systems 
were  restricted  to  known  stable  systems,  and  with  driving 
functions  far  from  resonance.   The  method  of  Flake  was  found 
to  be  so  impractical  that  it  is  only  mentioned  here  for 
reference. 


II.   METHODS 

The  various  methods  of  obtaining  approximate  solutions  to 
nonlinear  differential  equations  were  investigated. 

A.   PIPES'  METHOD 

Pipes  has  postulated  a  general  nonlinear  differential 
equation  of  the  form 

where  t  is  the  independent  variable,  x  is  the  dependent  vari- 
able, k  is  a  constant,  0 (t)  is  a  driving  function,  and  the  a. 
are  functions  of  the  differential  operator  D,  where 

D  -  d/dt  (2.1.2) 

Consider  the  special  case  of  equation  (2.1.1)  where 

and  the  a.  are  constant  coefficients  for  i  =  2,  3,  ...,  n. 
Equation  (2.1.1)  becomes 

1  cTt   Z<U*         dW        n       n  ,   .  .        (2.1.4) 

or  using  the  operator  notation 


CW+kD+^0%  .  ..^bnDn)*  +■  °^a  +  •  *  • 


(2.1.5) 


Following  Pipes,  assume  a  series  solution  of  the  form 

X(i)  =  A^W  +  AcokSAmk3*  . . .  (2.1.6) 


By  substituting  equation  (2.1.6)  into  (2.1.1)  and  equating 
coefficients  of  equal  powers  of  k,  the  various  A. (t)  terms 
may  be  found .   They  are 

A,W-  tftft/a,  (2.1.7) 

/\2ttW    -  a^AV^/a,  (2.1.8) 

A3(o  =  ~  Wot^*M>+  a5^i^l  (2.i.9) 

A4W  =  -  0/a)[aAA>^V^)^c^5^^A^UaA>]  (2.1.10) 


+  4ol4^i^K1Wa^K^3 


(2.1.11) 


Pipes  (Ref.  1)  lists  additional  terms.   However,  since  a 
second  or  third  approximation  to  the  solution  will  generally 
be  all  that  is  desired,  only  the  first  five  terms  are  given 
here. 

A  first  order  approximation  can  be  found  by  solving  the 
linear  differential  equation  (2.1.7).   Since  (2.1.7)  is  a 
linear  differential  equation,  it  can  be  solved  by  the  usual 
techniques  of  integration  or  Laplace  transforms,  applying 
initial  conditions.   The  remaining  A. (t)  are  solved  in  like 
manner,  except  the  initial  conditions  are  zero. 

B.   BRADY-BAYCURA  METHOD 

Brady  and  Baycura  have  shown  that  for  exponential 
functions  and  for  small  values  of  time  the  following  approxi- 
mate expression  holds 


Since  the  Laplace  transform  of  0(t)    is  known,  all  trans- 
forms form  the  terms  of  equation  (2.1.3)  are  known  a^d  thus 
(2.1.3)  can  be  solved  for  X(s). 

C.   TOU-DOETSCH-PIPES-NOWACKI  METHOD 

This  method  probably  originated  with  Pipes,  but  was 
arrived  at  independently  by  the  others. 

Assume  the  nonlinear  system  differential  equation  of  the 
form 

tantf  *  S-^  ...U.^  ...  -vex,  \y   Ri^.y.-O  =  q(0      (2.3.1) 

Taking  the  Laplace  transforms  of  both  sides  of  equation 
(2.3.1)  yields 

V^  Iavs^Qun(\K^^r-^^         <2-3-2) 

where  Q(s)  is  the  set  of  initial  conditions. 

In  the  first  approximation  assume  the  nonlinear  term 

i&VvV-il"0'  then 

VW.   Q»>  "  §  S  ■  (2.3.3) 

±  ^K 

The  inverse  Laplace  transform  of  Y,  (s)  yields  the  time 
domain  solution  of  the  first  approximation. 

The  nonlinear  term  is  included  in  the  second  approximation. 

X«-    ^^^,v---'J  (234) 
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where  f 


(t,  y,  /  Y-,  /  ...)  is  a  function  of  the  first  approxi- 


mation of  y, (t) . 

Successive  approximations  are  found  by  continuing  the 
same  procedure. 


III.   APPLICATIONS 

Nonlinear  systems  investigated  were  limited  to  second 
order  with  initial  conditions  and  extended  to  systems  with 
forcing  functions. 

A.   NONLINEAR  SPRING 

Consider  a  simple  spring-mass  system  described  by 

m  cN    +  F<w  =  o  (3.1.D 

The  system  is  undamped,  with  an  initial  displacement  of 
x(0)  =  x  ,  and  initial  velocity  set  to  zero.   F(x)  is  the 
spring  force,  and  for  this  example,  let  its  characteristic 
be 

Fo^>  =  kx  -V  bxZ  (3.1.2) 

Letting  all  constants  be  unity,  the  system  equation  is 
nonlinear 

cLA  4.  x  -V  XZ  -  6  (3.1.3) 

cH1- 


(3.1.4) 


1.   Solution  by  Pipes  Reversion  Method 

Using  the  method  outlined  in  Chapter  II,  Section  A, 
equation  (3.1.4)  can  be  rewritten  as 

10 


(3.1.5) 


The  a.  coefficients  are  seen  to  be 

1 

&x  K     \  \  (3.1.6) 

The  individual  A.  coefficients  are  then  computed  using  the 
formulas  (2.1.8),  (2.1.9)  and  (2.1.10): 

A,<*\«*<W/ft,  (2.!.,) 


Taking  the  Laplace  transform  of  equation  (3.1.3)  and 
impressing  initial  conditions  yields: 

A,<$=  iTTi  (3.1.8) 


This  well  known  transform  gives 

AtCO*  x^co^t  (3.1.9) 

As  expected,  the  system  is  oscillatory. 

The  next  term  A?(t)  is  found  from  equation  (2.1.8) 

A*W*  ""^   '    VH  (3.1.10) 

A2(t)  and  all  subsequent  A±(t)  are  determined  by  the  con- 
dition that  they  have  zero  initial  conditions. 
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Using  the  identity 

COS7 -I  -  ^0*0**0  (3.1.11) 

and  taking  Laplace  transforms  of  both  sides  of  equation 
(3.1.10) 

Az^=  'xUcFS      0b^.)(5^4T~  J  (3.1.12) 

Both  terms  can  be  expanded  into  partial  fractions : 

(3.1.13) 


Taking  inverse  Laplace  transforms : 

\(&  =  ~~P  "  3°^  _3  CosZi  J  (3.1.14) 

A~(t)  is  found  in  a  similar  manner,  and  the  solution  is 
given  by 

Xtfc,:  ^wV+A^\f  jv  ^jWk\  .  .  ,  (2.1.6) 

"*M  ^isvo-l-  WilCod-vlCt»b2i^CjDS3-y^..  (3.1.15) 
3  L  *■  ^6      3       it     -J    v 

In  equation  (3.1.15)  it  is  observed  that  the  first 
secular   term  occurs  in  the  third  approximation. 


2 
1Secular  terms  are  those  involving  t  sin  t,  t  cos  t,  t   sin  t, 

• . . /  etc . 
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2.   Solution  by  Brady-Baycura  Transform  Method 

Taking  the  Laplace  transform  of  equation  (3.1.3)  and 
applying  initial  conditions  gives 

5?XcO  -*X*  vXteb  \   5>X\0  ^=0  (3.1.16) 

In  the  usual  manner,  equation  (3.1.16)  is  solved  as  a 
quadratic.   After  rearranging  terms 


Neglecting  the  negative  radical  as  physically  impossible,  the 
positive  radical  can  be  expanded  as  a  binomial 

**'    ^^    *?**  (3.1.3.8) 


The  first  term  is  just  x   cos  t;  however^  L(sVh* 

is  not  one  of  the  common  transforms.   Reference  10  lists  the 
following  transform: 


71— 


^rM,,^ 


(>He£)nJ     2n"y,,rC^^n"3A  (3.1.19) 


where  Jn-3/2(w  t)  is  the  half  order  Bessel  function. 

Further  it  is  known  from  the  theory  of  Laplace  trans- 
forms that 
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Z  ly^iM  <--« +z  + 


G-Wa    ,  «-».»,«,'•■  ,3.1.20) 


(k) 
Since  fv  ' (0)  -  0  for  k  -  0,  1,  2,  .  ..,  n-1,  the  second  term 

can  be  found : 

letting  <36=  \ 


V-WV-^-"*) 


(3.1.21) 


and  after  carrying  out  the  differentiations  of  equation 
(3.1.20),  the  second  term  of  equation  (3.1.18)  becomes 


"  **  Istswt  4  "t2Coi>\;  I 


(3.1.22) 


Repeated  application  of  the  method  will  yield  the 
third  term. 

The  approximate  solution  to  three  terms  is  the  truncated 
series 

^[|(^t-^s\>|t^i  ,  £«lf*  ■  '  '   (3.1.23, 
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3.   Solution  by  Tou-Doetsch-Pipes-Nowacki  Method 

Putting  equation  (3.1.3)  in  the  form  of  (2.3.1) 


Q  Y  \   ax  -V  K7"  =  Q 


(3.1.24) 


where 


=  C> 


J 


(3.1.25) 


Taking  Laplace  transforms  of  both  sides 


(3.1.26) 


Xc^%Q-^^^^l  = 


o 


(3.1.27) 


Solving  for  X,  (s)  as  in  equation  (2.3.2) 


X.(3 


*^ 


SS\ 


(3.1.28) 


Then  taking  inverse  Laplace  transforms  yields  the  first 
approximate  solution: 


\(^=  Xbcob-t 


(3.1.29) 
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The  second  approximation  is  found  by  applying  equation 
(2.3.3) 


xc»-  ^^-xfrw,,...)! 


Z^   aKs  (3.1.30) 


k= 


& 


*&-     ^  ~  £rlZ[*C~H]l  (3.X.31, 


Using  the  identity  of  equation  (3.1.11),  equation 
(3.1.31)  becomes 

X<A  -  4^-  -  -*»   f-L  4  -1—1  (3  1  32) 


Taking  the  inverse  Laplace  transform  gives  the  second  approxi- 
mation as 


*fc4V>VfiA->£j  \-  *.c*a  -1  c^zVj 


(3.1.33) 


Repeating  the  method  gives  x3(t)  and  the  result  is  the  same 
as  equation  (3.1.15). 

While  this  method  appears  to  be  different  than  Pipes 
method  of  Chapter  II,  Section  A,  consideration  of  the  manipu- 
lation of  the  solution  shows  the  two  methods  to  be  equivalent 
Since  the  method  of  Chapter  II,  Section  A  is  less  cumbersom, 

it  will  be  used  hereafter. 

16 


4.   Comparing  the  Methods 

Figures  1  and  2  are  phase  portraits  of  the  solution 
as  given  by  the  IBM  360  computer.   Runge-Kutta  integration 
technique  was  used  to  obtain  the  solution.   From  these  plots, 
it  is  clear  that  the  system  is  oscillatory  for  x  ^0.5  units 
and  unstable  otherwise. 

Figures  3,  4  and  5  compare  the  solutions  by  Pipes 
method  and  the  Brady-Baycura  technique  for  various  initial 
conditions. 

The  graphs  demonstrate  the  relative  accuracy  of  the 
approximations.   The  Brady-Baycura  technic  gives  a  better 
approximation  in  each  instance;  however,  Pipes  method  is  as 
accurate  for  small  initial  conditions.   The  difference  is 
caused  by  the  location  of  the  secular  terms  in  the  solution. 
For  both  methods,  however,  the  approximations  are  of  little 
value  for  times  greater  than  one  unit. 

No  attempt  has  been  made  to  suppress  the  secular  terms 
which  cause  the  approximate  solution  to  diverge.   It  is  gen- 
erally known  that  subharmonic  resonance  terms  (submultiples 
of  the  natural  frequency)  exist  in  these  nonlinear  systems; 
they  do  not  occur  in  this  solution.   The  higher  frequency  jump 
resonance  terms  do  occur. 

B.   NONLINEAR  SPRING  WITH  DAMPENING 

The  system  of  example  3.1  was  modified  to  include  frictional 
dampening.   The  resulting  system  equation  takes  on  a  velocity 
term.   Letting  all  coefficients  be  unity,  except  the  velocity 
term,  the  system  equation  becomes 
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FIGURE  1.   Phase  portrait  of  x  +  x  +  x 
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FIGURE  2.   Phase  portrait  of  x  +  x  +  x 
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TABLE  I 

x   = 
o 

1 

Time 

Runqe-Kutta 

Brady 

Pipes 

0.000 

1.00000 

1.00000 

1.00000 

0.050 

0.99750 

0.99750 

0.99854 

0.100 

0.99002 

0.99002 

0.99418 

0.150 

0.97762 

0.97762 

0.98697 

0.200 

0.96039 

0.96038 

0.96424 

0.500 

0.76476 

0.76469 

0.86483 

1.000 

0.20173 

0.22176 

0.56088 

1.500 

-0.43604 

-0.23227 
0.50 

0.24276 

0.000 

0.50000 

0.50000 

0.50000 

0.050 

0.49906 

0.49906 

0.49919 

0.100 

0.49626 

0.49626 

0.49678 

0.300 

0.46675 

0.46674 

0.47137 

0.500 

0.41002 

0.41002 

0.42252 

1.000 

0.17967 

0.18245 

0.22407 

1.5000 

-0.10512 

xo  = 

-0.075097 
0.25 

-0.02327 

0.000 

0.25000 

0.25000 

0.25000 

0.050 

0.24961 

0.24961 

0.24963 

0.1000 

0.24844 

0.24844 

0.24850 

0.500 

0.21213 

0.21213 

0.21369 

1.000 

0.11163 

0.11214 

0.11715 
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cHl   ^  cH  (3.2.1) 


(3.2.2) 


The  investigation  included  solving  the  system  equation  by 
Pipes  and  Brady-Baycura  techniques.   The  final  solution  is 
found  by  letting  the  coefficient  C  approach  zero  in  both 
methods. 

1.   Solution  by  Pipes  Method 

Using  operator  notation,  the  a.'s  are 

a  *  if  -*  et>  -v  \ 


(3.2.3) 


By  letting  C  =  2  A  ,  the  algebra  is  simplified  and  the  same 
result  will  obtain  by  letting  2 fi    go  to  zero. 

Proceeding  with  the  solution  A,  =  0  /   a,  and  upon 
substituting  and  rearranging 

CD%CtuO\<^S  (3.2.4) 

Taking  Laplace  transforms  and  impressing  initial  conditions, 
equation  (3.2.4)  becomes 
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C^tStOM^H  (3.2.5) 


or 


J\   (^  =   x^mQ 

>       (&fl-Tfi&)(**fl+ygKT)  (3.2.6) 


Reference  11,  formula  119,  gives  the  inverse  trans- 
form of  equation  (3.2.6)  and 


(3.2.7) 


letting    2  A    go   to   zero 

AC-0-    K.^v>^3  -VexpC-jO]  (3-2-8) 


'    .  (3.2.8a) 


The  solution  for  A_ ,  after  a  modicum  of  algebra 
becomes 

\M-  -^^\-%<*^-^C*«V]  (3.2.9) 


and  substituting  A^(t)  into  equation  (2.1.5),  Ao(t)  becomes 


3   T?*cT>xl 


(3.2.10) 
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or 


Ai^  tt+MI»»L*  J      (3-2-11} 


Taking  Laplace  transforms,  equation  (3.2.11)  becomes 

A  ®  -.    <  \ 5 +  5^ 

(3.2.12) 

v 


Va*  ^s 


(^-v^%2/6^+  0    ^^(^  *  Z^if  & 


All  inverse  transforms  of  the  terms  of  equation  (3.2.12), 
except  the  second,  can  be  found  in  Ref.  10.   The  second  term 
is  evaluated  by  convolution  as  follows: 
Let 


s  -*-V  (3.2.13) 


^°  ~       sMy^+l  (3.2.14) 


then   z1 (t)    ~   cos   t   and    z?(t)    -   sin  t  when   2/6   goes   to   zero. 
Convolving    z,    and    z~ 

H^ta  *  Z^Ctt  -  J   W  2-z  (*-■*)  oft. 

c     t.   Sxo-t  (3.2.15) 


The  approximate  solution  to  three  terms  by  Pipes 
technique  is 
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K  to  X     X.CD^i  -  X 


■I- 


"^^  ~4  i'voi-^,  v\M  CWsl~i-^S2t-^-CjA3i)-V  ••  .    (3.2.16) 
r6  \  +   W  4-8.       «       »»=»_] 


It  is  observed  that  the  first  secular  term  occurs  in 
the  third  approximation. 

2.   Solution  by  Brady-Baycura  Method 

The  solution  proceeds  by  taking  the  Laplace  trans- 
form of  equation  (3.2.1),  and  impressing  the  initial  con- 
ditions.  The  nonlinear  Laplace  transform  is  used  for  the 

2 

x  (t)  term. 

^X^UC^^+CX^  =  **(Mc}  (3.2.17) 

By  reversion  of  series  technique,  the  left  side  of 
equation  (3.2.17)  can  be  made  into  a  series  by  using  the 
formulas  of  Ref.  9. 

Let  equation  (3.2.17)  have  the  form 

*fcCs»tcV  ^Xc^V^X^-v  ,  %  .  =  u  (3.2.18) 

where 

I     ^       f  (3.2.19) 
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then  by  reversion  of  the  series 

X^V-  ^V^^vP>^\.  ,  ,  ,  (3.2.20) 


where 


-J    to  (3.2.21) 


B,  =   b.      4X^M)  (3.2.22) 


&.«- 


L 


b?        C?^^+Iy  (3.2.23) 


b3,  ^(*SAO-    j^ 


or  letting  C  -  2/2> 


(3.2.24) 


^     2*^(^0*  ^   |t(  (3.2.25) 

(^■•VZ//SSArO£' 


The  first  term  of  equation  (3.2.25)  is  the  same  as 
equation  (3.2.6),  and  is  just  x   cos  t,  as  in  the  Pipes 
solution. 

Evaluation  of  the  second  term  must  be  found  by  con- 
volution : 

Let 
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■2:  (s)  - 

g-+Zfi*t  \  (3.2.26) 


-^  -  -^liL____  (3.2.27) 


??-*  ^^S  -v  i 


Taking  inverse  transforms,  and  letting  2a  go  to  zero, 

"i^55    ^^  •  (3.2.28) 

^TCt^  =•  eo<bt  (3.2.29) 


Then  the  inverse  of  the  second  term  is  the  con- 
volution 

2ct^  -  ^,(.-0  ^-2:^  £  ^C-O  (3.2.30) 

Again  after  a  modicum  of  calculus,  algebra  and  trigonometric 
substitutions,  equation  (3.2.30)  reduces  to  just 


2?<«  =    "t  \  V^^  v~^o-t  j 


(3.2.31) 


and  to  two  terms,  the  Brady-Baycura  technique  yields 

XCO  -x.  <D^os^  -  *M  5  cdcA.  V-ti^Vl  >  ,  .  .       (3.2.32) 
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The  technique  has  lost  its  desirability  even  at  this 
point,  since  the  inverse  transforms  are  so  difficult  and  time 
consuming  to  evaluate. 

3.   Comparing  the  Solutions 

The  phase  portraits  given  in  Figures  6,  7a  and  7b 
show  the  apparent  warping  of  the  phase  plane  due  to  the 
inclusion  of  the  velocity  term.   The  instability  of  the 
system  is  seen  to  depend  on  both  the  value  of  C  and  x  . 

The  Runge-Kutta  solution  is  compared  with  the  approxi- 
mate solutions  of  equations  (3.2.16)  and  (3.2.32)  in  Figures 
8,  9  and  10.   Table  II  lists  selected  values  of  the  Runge- 
Kutta  solution  and  the  approximations. 

It  is  clear  from  the  table  that  neither  approximations 
are  good  even  at  0.5  time  units.   It  is  further  obvious  that 
for  values  of  the  coefficient  CLL1.0,    the  system  equation 
(3.2.1)  will  approach  equation  (3.1.3).   Under  the  condition 
of  C  small,  the  techniques  might  have  some  value  in  approxi- 
mating the  transient  state  for  small  times  (milliseconds) . 

C.   NONLINEAR  SPRING  WITH  CUBIC  TERM 

Consider  again  equation  (3.1.1),  but  let 

Fto  -  t*  -V  W.3,  (3.3.1) 

Also  let  all  coefficients  be  unity  and  the  system  have  the 
same  initial  conditions  as  before.   The  nonlinear  expression 
describing  the  system  is 


£jL   \  X  - 


a-tv  (3.3.2) 
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FIGURE  7a.   Phase  portrait  of  x  +  0.5x  +  x  +  x     0 ,  xq 
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32 


FIGURE  7b, 


Phase  portrait  of  x  +  0.5x  +  x  +  x 


2  _ 


0,  x 


1.2 
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TABLE  II 

X   = 

o 

1.0 

Time 

Runqe-Kutta 

Brady 

Pipes 

0.000 

1.00000 

1.00000 

1.00000 

0.050 

0.99754 

0.99719 

0.99667 

0.100 

0.99035 

0.98877 

0.98669 

0.150 

0.97870 

0.97478 

0.97014 

0.200 

0.96290 

X 

o 

0.95530 
0.50 

0.94711 

0.000 

0.50000 

0.50000 

0.50000 

0.050 

0.49908 

0.49898 

0.49896 

0.100 

0.49638 

0.49594 

0.49584 

0.150 

0.49200 

0.49089 

0.49066 

0.200 

0.48604 

X   = 

o 

0.48384 
0.25 

0.48344 

0.000 

0.25000 

0.25000 

0.25000 

0.050 

0.24962 

0.24959 

0.24960 

0.100 

0.28849 

0.24836 

0.24839 

0.150 

0.24666 

0.24632 

0.24638 

0.200 

0.24417 

0.24347 

0.24357 
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(3.3.3) 


1.   Solution  by  Pipes  Technique 

Using  operator  notation,  the  a.  coefficients  are 
identified  as 


a ,  =  IN  \ 

CX^--    D 

^»-  N 

V  -  » 

QSC«=  fe 

(3.3.4) 


The  A.  (t)  coefficients  are  found  next.   From  equation  (2.1.7) 


or 


(X,%OA,W^  _       13.3.5, 

Taking  Laplace  transforms  and  impressing  the  initial 
conditions  on  A, (t)  the  solution  is  immediately  seen  to  be 
the  same  as  that  of  the  quadratic  case 

AJ^X.CdA  (3.3.6) 


Substitution  of  a2  into  equation  (2.1.4),  gives  coefficient 
A_  as  zero.   A~  is  found  from  equation  (2.1.4): 
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A5«  ■  ■■$«•&*& +  ^>\] 


(2.1.8) 


\to-4nb  ♦*****■ 


(3.3.7) 


3 

and  after  substituting  for  cos   t,  equation  (3.3.7)  becomes 

A^^L^r^l  (3-3-8) 


Taking  Laplace  transforms  of  both  sides 


A^  =    *?  \   _^    ^       s       |  (3.3.9) 


*3 


4-  L  ^4  OP-    CiS  ^  (j^  -v  o 


The  inverse  transforms  for  the  factors  of  equation  (3.3.9) 
are  found  in  reference  11,  and  the  solution  of  A^(t)  is 


J     "S"  A       ^V       J 


Since  all  even  a.  and  thus  even  A.  are  zero,  A.  is 

l  i  4 

zero.   The  coefficient  Ac  ,  after  substitution,  reduces  to 


AEto=-^[^KM^] 


(3.3.11) 
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Substituting  for  the  coefficients  and  manipulating  trigono- 
metric identities, 

32.^0 


4-         <V 


(3.3.12) 


and  taking  Laplace  transforms 


~^r  L^-^  &+0        ^o5        v?x*o* 


-  'A*  VA^>      1        (3.3.13) 


(P^C^vO    £ 


The  last  three  terms  have  readily  available  inverse  trans- 
forms; however,  the  first  and  second  terms  require  special 
techniques. 

First  term  is  evaluated  by  convolution: 

Let 


*.w«£  [  sfe-]*  c&^ 


(3.3.13a) 


^C-0  -  \.  L  k>  MiS<o  J  : 


^<K^^J  Z^\~  (3.3.13b) 


then 
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(3.3.13c) 


-*>■)_  °\(o  *A  (3.3.13d) 


The   evaluation   of    the    second    term    is    found   in  Appendix  A: 

•t'bS?]  "-    H^"^  1  ,3.3.13.) 

The  remaining  terms  of  (3.3.13)  are  found  in  reference  10. 
The  result  for  A,,  (t)  becomes 

^    ^<+        **      _\     (3.3.14) 


The  approximate  solution  by  Pipes  method  is  the  truncated 
series 

-^7|_S^)  J  (3.3.15) 


+  *(&v^-.3_=>%<*3*-"^fc>c*s-t 
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It  is  seen  from  equation  (3.3.15)  that  the  first 
secular  term  occurs  in  the  second  approximation. 
2.   Solution  by  Brady-Baycura  Method 

Upon  taking  Laplace  transforms  and  using  the  non- 
linear expression  of  equation  (2.2.2)  and  impressing  the 
initial  conditions,  equation  (3.3.2)  becomes 

2fxVa*tc-4r}Xo\  =  *K.  0.3.16) 

Considering  equation  (3.3.16)  to  be  a  series,  the  method  of 
reversion  outlined  in  Chapter  III,  Section  B,  paragraph  2 
can  be  used  here  also  to  obtain  all  the  terms: 


l-V  t>,      ^a-\  (3.3.17a) 


Bt  -  O  (3.3.17b) 

E^-ifC^VO--^-  ~~T  (3.3.17c) 

The    solution   in   terms    of   X(s)    is 

AC3^        a?-+i  (j^'V*  f^v.V  ~    "  '  (3.3.18) 


^0 


The  inverse  transform  of  the  first  term  is  just  xq  cos  t. 
The  second  term  must  be  solved  using  the  technique  outlined 
in  Brady  (Ref.  2)  and  in  the  example  of  Chapter  III,  Section 
A,  paragraph  2.   The  method  is  straightforward,  giving  the 
second  term  as 
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-  **>  \~\£-t£»v\\  -"^?5>v^^  "V  ^  -Ccosd. 


<=V-§>  L. 


(3.3.19) 


The  second  term  (3.3.19)  yields  the  higher  order  secular 

2  3 

terms  t   cos  t  and  t   sin  t  which  will  be  small  contributions 

for  t  LI. 

It  is  also  known  from  equation  (3.1.19)  that 

» 

■V  L  ^  ,y  j       Z'K  p^  _  (  (3. 3. 20) 


By  observing  the  result  in  second  term,  one  sees  that 
even  higher  order  (t  )  secular  terms  will  result  from  the 
third  term  of  equation  (3.3.18),  and  further  that  the 

coefficient  for  this  third  term  is  small  (less  than 

-4 
.2  x  10   ) ;  therefore,  little  is  gained  by  solving  equation 

(3.3.20) . 

The  approximate  solution  by  the  Brady-Baycura 
technique  is 

X<£*  X.cx*>1:  -  ^[^la-iU^U^,  (3.3.21) 


3.   Comparison  of  Solutions 

Figures  11  and  12  show  the  phase  plane  solution  for 

equation  (3.3.2).   Here  the  solution  is  seen  to  be  stable  for 

all  x  .   The  system  always  remains  in  a  limit  cycle, 
o 
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It  is  also  observed  from  Figures  13,  14  and  15  that 
for  larger  initial  conditions,  the  Brady-Baycura  method  gives 
a  better  approximation  to  the  solution  for  early  time  inter- 
vals.  For  the  case  of  x   small,  Pipe's  method  gives  better 
results.   This  is  caused  by  the  location  and  number  of 
secular  terms  appearing  in  the  solution,  which  causes  it 
to  diverge. 

Selected  values  for  the  graphs  of  Figures  13,  14  and 
15  are  given  in  Table  III. 

D.   NONLINEAR  SPRING  WITH  DRIVING  FUNCTION 

Consider  the  system  of  Chapter  III,  Section  A,  but  now 
as  a  forced  system.   The  system  nonlinear  differential 
equation  is 

cRv  (3.4.1) 

where 

v^=W^W-  <3-4-2> 

Letting  all  coefficients  be  unity,  equation  (3.4.1)  becomes 
°^-^^'  -  "aO^iV  (3.4.3) 


(3.4.4) 
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FIGURE    11.       Phase   portrait   ofx+x+x  0,    x  1.0 
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FIGURE  12.   Phase  portrait  of  x  +  x  +  x 


3 
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FIGURE  13.   Time  s 
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FIGURE    15.       Time    solution   of   x   +   x   +  x  0 ,    xQ         0.25 
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TABLE 

III 

x   =  ] 
o 

..0 

Time 

Runqe-Kutta 

Pipes 

Brady 

0.000 

1.00000 

1.00000 

1.00000 

0.050 

0.99750 

0.99768 

0.99750 

0.100 

0.99003 

0.99073 

0.99002 

0.150 

0.97767 

0.97920 

0.97761 

0.250 

0.93877 

0.94277 

0.93831 

0.500 

0.76880 

0.78095 

0.76278 

1.000 

0.23370 

x   -  0. 
o 

0.25179 
5 

0.19357 

0.000 

0.50000 

0.50000 

0.50000 

0.050 

0.49922 

0.49922 

0.49922 

0.100 

0.49688 

0.49690 

0.49688 

0.250 

0.48064 

0.48077 

0.48063 

0.500 

0.42461 

0.42498 

0.42444 

1.000 

0.22682 

0.22707 

0.22681 

x   =  0. 
o 

,25 

0.000 

0.25000 

0.25000 

0.25000 

0.050 

0.24967 

0.24967 

0.24967 

0.300 

0.23815 

0.23816 

0.23815 

0.750 

0.17928 

0.17930 

0.17931 

1.000 

0.12946 

0.12947 

0.12966 

1.500 

0.0089567 

0.0089509 

0.011029 
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1.   Solution  by  Pipes  Method 

The  a.  coefficients  are  those  of  equation  (3.1.6), 
except  here  it  is  seen  that  0 (t)  =  B  cos t^  t. 


The  solution  proceeds  as  before 


^>  <3-4'5> 


Taking  Laplace  transforms  and  impressing  initial  conditions, 


AbV-  ■  ?^    +  — ^ (3.4.6) 


The  terms  of  A, (s)  have  transforms  in  Ref .  10 


(3.4.7) 


W,CD<d.  -W.Cj^cM;  (3.4.8) 


where 


Wx=  ^   *  _Jk 


^"-\ 


(3.4.8a) 


The  solutions  to  A2(t)  and  A3(t)  are  straighf orward 
using  the  method  outlined  in  Chapter  II,  Section  A.   Thus 

*l^*t^-^<^«<2tfT[  (3.4.9) 
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where 


\u-      ^ 


-  ^ 

X*^-) 


(3.4.9a) 


.7 


The  solution  for  A_(t)  is  given  in  Appendix  B. 
The  complete  solution  by  Pipes  method  is  given  by 


XIO*  A.^-0  v  kl-0  v  ^l*U  ,  .  » 


(3.4.10) 


2.   Solution  by  Brady-Baycura  Technique 

Equation  (3.4.3)  can  be  transformed  as  in  Chapter 

III,  Section  A,  paragraph  2  using  the  nonlinear  transform 

2 

for  x  (t) .   After  impressing  the  initial  conditions,  the 

transformed  system  equation  becomes 


5>X  C^  \  C^v^X  e^  =  *,>  -v  1^  ^ 


^-v^ 


(3.4.11) 


Letting  equation  (3.4.11)  be  a  series,  then  by  the  method  of 
algebraic  series  reversion 


-\ 


(3.4.12) 


where  the  A.  coefficients  are 

l 
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Al. 


^>  +  \ 


5. 


(S--V  s$ 


(3.4.13) 


,i_ 


Thus 


&\\ 


iS^ 


C^+»p\_     ^^j 


(3.4.14) 


The  first  term  of  the  series  in  equation  (3.4.14) 
is  easily  evaluated  by  performing  the  algebra  indicated  and 
using  the  transform  pairs  of  Ref.  10.   Let  the  first  term 


be  P. (s) ,  then 


F(^- 


*»^        "^« 


S^iri      ^0V??-*«*t^  (3.4.15) 


and 


\&  "  *+**&    v  -S-W\  -  <Los^4, 1  (3.4.16) 


-  b^cjosi  -V,oe>st5d. 


(3.4.17) 


where 
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V 

w^ 


*o   t 


-& 


U^-v 


"fe 


«*-! 


(3.1.18) 


The  result  is  not  unlike  the  first  approximation 
obtained  in  the  previous  method  by  Pipes. 
Let  the  second  term  be  F?(s),  then 


^ 


*&' 


<*^ 


^ 


(3.4.19) 


r   ■«- 


*6£>- 


Z&X^" 


— ■»"*-  T>      — i 


L  (^  0*    <*\  O"5  OHaf)   Cs?-+i^(sH^0 


(3.4.20) 


The  terms  of  F2(s)  are  all  uncommon  transforms.   The 
technique  of  finding  the  inverse  transforms  is  as  follows: 
Since 


^  IjP^?  J   »  ^"k-  Vte^) 


(3.4.21) 


As  shown  in  Appendix  A,  the  derivative  method  of  finding 
inverse  transforms  is  extended  to  equation  (3.4.21). 


— (U>^- W)]  -  -W\^^w^] 


(3.1.21a) 


and  the  inverse  transform  of  the  first  term  of  equation 

(3.4.20)  is 
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J^Ta-taiA*  -v -\>ex>sVj 


(3.4.22) 


The  second  term  of  F,,(s)  is  evaluated  by  convolution.   Except 
for  a  constant,  it  is  the  first  term  of  f~(t)  convolved  with 


Cbi 


^ 


or 


*fcxj)fel?ts*nt  V-^CLDd:    4  W^*\t 


(3.4.23) 


The   expression    (3.4.23)    can  be   evaluated    as 


W«v  \_ 


•a: 


•\V 


.    (oi.       ^vnfeuy^    -  J5L_CwJt«JV\)t 


O0-N 


(<A>-V 


(3.4.24) 


where 


c 

*.v4= 

CiV  tf- 

tw- 

u 

^>-o" 

Vl.« 

<*> 

\(o 

UV*1 

w.,- 

z 

4 


<^«-03      c^vo* 


4 


4 


ce^o"     csf^»y    <h-o3 

Vd  _        4 .^         4 


t*V 


c^-v  Qu^D"        l!A,-^ 


^v\ 


(3.4.24a) 
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Except  for  a  constant,  the  last  term  of  F  (s)  is 
simply  the  term  of  (3.4.24)  convolved  again  with  1_  ^wnuN-V 
The  result  is  found  in  Appendix  C. 

Clearly  the  method  is  too  cumbersome  to  carry  beyond 
this  second  approximation.   The  approximate  solution  by  the 
Brady-Baycura  technique  is 

*«**,«*V&*...  (3.4.26) 


3.   Comparing  the  Solutions 

The  time  plots  of  Figures  16,  17,  18  and  19  exhibit 

the  relative  accuracy  of  the  approximations  to  the  Runge- 

Kutta  solution,  for  various  combinations  of  x  ,  B  and  C?)  . 

o  o 

From  Table  IV  it  is  apparent  that  the  approximations 
are  good  only  for  the  first  initial  time  units.   Pipes  method 
gives  only  a  slightly  better  approximation  than  the  Brady- 
Baycura  technique.   For  the  case  uS  '     =   5.0,  it  is  observed 
that  the  approximations  by  both  methods  improve.   This  is 
in  keeping  with  the  original  restriction  of  the  system  to 
forcing  functions  far  from  resonance. 

E.   NONLINEAR  SPRING  WITH  CUBIC  TERM  AND  FORCING  FUNCTION 
Assume  the  describing  equation  to  be 

^  (3.5.1) 

with  initial  conditions 
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Time    solution   of   x   +  x   +   x      -0.5    cos    1.5t, 


x      =    0.5 
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18.       Solution   to  x   +  x   +  x      -    0.5    cos    1.5t,    x  0.25 
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Time    solution   of   x   +   x   +   x      "0.25    cos    1.5t, 


x      =    0.5 
o 
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Time 


0. 

,000 

0. 

,050 

0. 

,100 

0. 

,250 

0. 

,500 

1. 

,000 

1. 

,500 

0.000 
0.050 
0.100 
0.250 
0.500 
1.000 
1.500 


0.000 
0.050 
0.100 
0.250 
0.500 
1.000 
1.500 


TABLE 

IV 

x   -  B 
o 

=  0.5, 

-  1.5 
o 

Runge-Ku 

itta 

Pipes 

Brady 

0.50000 

0.50000 

0.50000 

0.49969 

0.49969 

0.50015 

0.49875 

0.49875 

0.50053 

0.49208 

0.49205 

0.50217 

0.46719 

0.46657 

0.49946 

0.35328 

0.34448 

0.41322 

0.13099 

0.09665 

0.14546 

x   =  B 
o 

=  0.5, 

=  5.0 
o 

0.50000 

0.50000 

0.50000 

0.49968 

0.49968  . 

0.49991 

0.49870 

0.49870 

0.49955 

0.49034 

0.49030 

0.49466 

0.44401 

0.44340 

0.45529 

0.17816 

0.17054 

0.17995 

•0.12313 

-0.14487 

-0.18418 

x   =  0. 
o 

25,  B  = 

;  0.5,     =  1.5 
o 

0.25000 

0.25000 

0.25000 

0.25023 

0.25023 

0.25041 

0.25093 

0.25093 

0.25161 

0.25563 

0.25563 

0.25960 

0.26987 

0.26979 

0.28343 

0.29115 

0.28972 

0.32447 

0.22963 

0.22239 

0.26635 

0.5,    B    =    0.25, 


1.5 


0.000 
0.050 
0.100 
0.150 
0.250 
0.500 
1.000 
1.500 


0.50000 
0.49937 
0.49750 
0.49438 
0.48444 
0.43861 
0.26682 
0.015343 


0.50000 
0.49938 
0.49750 
0.49438 
0.48441 
0.43801 
0.25904 
-0.01124 


0.50000 
0.49972 
0.49881 
0.49722 
0.49164 
0.46008 
0.29380 
•0.01725 
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XCpV  X, 


(3.5.1a) 


1. 


a.  '  s. 


Solution  by  Pipes  Method 

The  solution  proceeds  in  the  usual  way  be  identifying 

t 


(3.5.2) 


Using  equation  (2.1.7),  the  first  approximation  A,,  is  found 
by  Laplace  transform  techniques. 


Am>>  e 


(3.5.3) 


Av^  - 


*o^ 


"&S 


■A-  \  «^"V\KsT-V^ 


(3.5.4) 


The  solution  is  the  same  as  for  the  quadratic  case 


A,C*V-  ^,C.O<sL-  Wo.CjOS.vx^t 


(3.5.5) 


where 


V 


R> 


(3.5.6) 


*£--) 
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The  coefficient  A  (t)  was  seen  previously  to  be  zero. 
Solving  for  A3(t)  is  straightforward  and  the  result  is 


\ i Cos>3>cob-i 


(3.5.7) 


where 


k*B   ^  V  >k,£ 


^ 


3 


^"-  %^K*  3/^x 


(3.5.8) 


The  solution  of  the  undriven  system  of  Chapter  III,  Section 
C  should  demonstrate  the  unfeasibility  of  Carrying  the 
approximation  beyond  A^(t).   Furthermore,  the  first  secular 
term  in  A_(t)  will  cause  the  solution  to  diverge. 
2.   Solution  by  Brady-Baycura  Technique 

Taking  Laplace  transforms  of  equation  (3.5.1)  and 
impressing  the  initial  conditions,  the  system  equation 
becomes 

^X5(>)  ^C^J^V'SXo-*— Jb (3.5.9) 
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Using  the  algebraic  reversion  of  series  technique  to  solve 
for  X(s),  let  the  left  side  of  equation  (3.5.9)  be  a  series 


k  =  -tav  t  xcsw  Wx*io  A  ^X^vK  ,  .  . 


(3.5.10) 


then 


^V^.^^N^V..  .  (3.5.u) 

=  ^,  k  ,  -^r-V- VsCV-Jl.  \v  .  .  .   (3.5.12) 


D 


Where  from  Ref.  9 


)  \ 


^B   V  "   5>V 


(3.5.13) 


1>.  -  -i 


=  D 


^  (3.5.14) 

T2>3-  _L_(zy-_LvO--_W_  =  -^  (3.5.15) 

Then  to  first  order  approximation 


*<**  V«i  (3.5.16) 


S>**    v      "^ 


^-v  \ 


c^^x*>o 


(3.5.17) 


and 

KC^  =  k^c^-t  -kj^u^-t  (3.5.18) 


64 


where 


(3.5.18a) 


As  expected,  the  result  is  equivalent  to  Pipes  method  for 
the  first  approximation. 

Since  B2  is  zero,  the  third  term,  F^(s),  of  the 
series  is  found 

(3.5.19) 


*uv  "$s* 


s>*5  V 


TS>! 


^^ 


*,V 


V&X  < 


s  sr 


3~&\>* 


($-*$         (£"?(&**)*      (^(^) 


+ 


— *3  s- 


(3.5.20) 


(Jv^)^Nx^ 


The  inverse  transform  of  the  first  term  as  given  in 
equation  (3.3.18). 


4-S 


ifffcanvt  *SlWt  -^s^-lj 


(3.5.21) 


Except  for  the  constants,  the  successive  terms  of  F^(s)  must 
be  evaluated  by  successive  convolution  of  —  sin^N   t  with 
the  right  side  of  equation  (3.5.21).   While  continuing  the 
evaluation  of  F-.(s)  will  lead  to  terms  of  cos  t,  modulation 
terms  and  higher  resonance  terms,  it  will  also  contribute 
numerous  secular  terms,  and  thus  cause  the  solution  to  di- 
verge more  rapidly.   Again  the  technique  has  lost  its  desir- 
ability by  providing  cumbersome  transforms. 
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Using  the  first  term  only  of  F  (s)  ,  the  solution  is 
given  by 


*t% 


-iV*t|*..,     <3'5-22) 


3.   Comparing  the  Solutions 

Figures  20,  21,  22  and  23  are  graphs  of  the  approxi- 
mate solutions  and  the  true  solution  for  various  combinations 

of  x  ,  B  anduA  .   Table  V  lists  selected  values  from  the 
o         o 

solutions. 

The  graphs  demonstrate  that  except  in  the  extreme 
case  of  Figure  20  where  x   =  B  and  &   =  1.5,  that  Pipes 
Method  provides  only  a  slightly  better  approximation  than 
the  Brady-Baycura  Method.   This  is  due  to  the  number  of 
terms  generated  by  the  Pipes  Method.   Figure  21  shows  that 
the  condition  of  the  system  operating  far  from  resonance, 
Pipes  Method  gives  a  better  solution  for  a  longer  period  of 
time.   Both  approximations  improve  as  x   and  B  become  small. 
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FIGURE    20.       Time    solution   of   x   +   x   +   x  0.5    cos    1.5t, 

x      =    0.5 
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FIGURE  21.   Time  so 


lution  of  x  +  x  +  x3  =  0.5  cos  5.0t, 


x   =  0.5 
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FIGURE  22.   Time  solution  to  x  +  x  +  x     0.5  cos  1.5t, 


x   =  0.25 
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FIGURE  23.   Time  solution  to  x  "•"  x 

x   =  0.5 
o 


<  +  x  +  x3  =  0.25  cos  1.5t, 
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TABLE  V 


Time 

0.000 
0.050 
0.100 
0.150 
0.250 
0.500 
1.000 
1.500 


x   -  B  - 

o 

:  0.5, 

=  .  1.5 

Runge-Kutta 

Pipes 

0.50000 

0.50000 

0.49984 

0.50017 

0.49937 

0.50068 

0.49857 

0.50150 

0.49594 

0.50386 

0.48209 

0.51016 

0.40483 

0.47304 

0.22075 

0.28681 

Brady 

0.50000 
0.49984 
0.49937 
0.49858 
0.49599 
0.48280 
0.41464 
0.26156 


-  B  -  0.50, 


=  5.0 


0.000 
0.050 
0.100 
0.150 
0.250 
0.500 
1.000 
1.500 


0.50000 
0.49984 
0.49932 
0.49834 
0.49420 
0.45886 
0.22743 
■0.04477 


0.50000 
0.49988 
0.49949 
0.49870 
0.49516 
0.46195 
0.23117 
-0.04985 


0.50000 
0.49984 
0.49933 
0.49835 
0.49425 
0.45941 
0.23216 
0.02377 


0.000 
0.050 
0.100 
0.250 
0.500 
1.000 
1.500 


x   ~  0.25,  B 


0.25000 
0.25029 
0.25116 
0.25709 
0.27572 
0.31407 
0.27665 


=  0.50, 


=  1.5 


0.25000 
0.25034 
0.25135 
0.25816 
0.27884 
0.31302 
0.25064 


0.25000 
0.25029 
0.25117 
0.25711 
0.27596 
0.31748 
0.29059 


0.000 
0.050 
0.100 
0.150 
0.250 
0.500 
1.000 
1.500 


x   ~  0.50,  B 
o 

0.50000 
0.49953 
0.49812 
0.49578 
0.48829 
0.45336 
0.31626 
0.09741 


0.25, 

0.50000 
0.49962 
0.49846 
0.49653 
0.49029 
0.45969 
0.32228 
0.07727 


=  1.5 


0.50000 
0.49953 
0.49813 
0.49579 
0.48831 
0.45362 
0.32072 
0.12178 
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IV.   CONCLUSIONS 

It  has  been  demonstrated  that  these  various  methods  of 
obtaining  approximate  solutions  to  nonlinear  systems  with 
initial  conditions  are  accurate  only  during  the  early  time 
intervals.   Further,  it  was  shown  that  this  degree  of  accuracy 
may  be  satisfactory  for  engineering  applications. 

It  has  also  been  shown  that  the  Brady-Baycura  technique 
of  nonlinear  transforms  usually  provides  a  more  accurate 
solution;  however,  the  ease  and  facility  of  Pipes  Method 
might  generally  make  this  a  more  desirable  method. 

The  method  of  handling  so-called  secular  terms  is  still 
unresolved.  .It  is  suggested  that  they  may  be  eliminated  by 
letting  the  frequency  also  be  a  series  such  that  the  con- 
stants of  the  series  are  chosen  to  eliminate  resonance  con- 
ditions and  thus  eliminate  the  secular  terms.   The  method 
would  be  similar  to  that  given  by  Pipes  in  Ref.  12. 

Both  methods  show  that  the  higher  order  harmonic  terms 
exist  in"  the  solution;  however,  the  subharmonic  terms  fail 
to  occur  in  the  first  few  terms  of  the  approximations. 

Since  the  problem  is  an  initial  value  problem,  it  appears 
that  the  solution  could  be  improved  by  stopping  the  problem 
solution  before  it  diverges  from  the  true  solution  (t   .5), 
then  restarting  with  a  new  set  of  initial  conditions.   The 
effects  of  x(o)  would  have  to  be  investigated.   The  value  of 
such  a  method  seems  dubious  since  Runge-Kutta  analysis  on  the 
digital  computer  seems  an  easier  method. 
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APPENDIX  A 


INVERSE  LAPLACE  TRANSFORM  OF 


(s2  +  l)3 


The  derivative  technique  is  used 
From  Ref.  10, 


(A-l) 


also 


A     _A 


4^ 


dsL^vOLJ  (sHO* 


(A- 2) 


or 


=  -±L_ 


(sH>)3  4ds[  (£>*y- 


(A-3) 


since 


X 


\    C.0 


-rMTFW  .  o.i.i.*...  , 


»  i   •> 


(A-4) 


where  n  is  the  n   derivative.   Thencombining  equations 


(A-4)    and    (A-3) 


rfcK 


»y 


4- 


■qi 


^ 


L 


CsHiT 


(A-5) 


iTsm^-  "ttjcvt 


(A-6) 
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APPENDIX  B 


EVALUATION  OF  A3  IN  FORCED  QUADRATIC  SYSTEM 


From  equation  (2.1.5)  A^  becomes 


(B-l) 


After  substituting  for  a..  ,  a?/  A   and  A   and  combining  terms, 

A3UV  _J H  c*>*A  -  W  W^ca^z\  -  W^  do^^-L  -  k  C3bb^ 


where 


(B-3) 
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Except  for  the  first  term  of  A^ ,  the  transforms  are  all  of 


the  form 


S 


(&S\}  UW) 


(B-4) 


For   which    the    inverse    transform    is  known 


-J \  <Ubs  V  -  to^tst 


(B-5) 


The  first  term  has  the  transform: 
which  has  the  inverse 

Using  expressions  (B-5)  and  (B-7) ,  A^  becomes 


(B-6) 


(B-7) 


A,W- 


-  W.W< 


W\h_  _  \z^  v  \jn  __  K\2_X^± 


r-   -v 


c^-o7--\     ov^-i     (3-^\y--\    (t-^-0->  •   ^^-\ 


+% 


(B-8) 


~  !^s  ■  c^s (to0-^  -b  _  W„   oxfepoK-tQ-b-  ^W  d&s>(co^  Ofc 


^G*o-0         4*oM>*v        ^(ft-0 


^-1 


-*.**] 
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APPENDIX  C 


SOLUTION  OF  LAST  TERM  OF  FJs)  BY  BPA.DY-BAYCURA  METHOD 

Let  the  expression  in  parenthesis  of  (3.4.24)  equal  Z, (t) 


and  —  sin     t  equal  Z~(t),  then 


1_ 

(o 


(C-l) 

thus 


^^L        J  (C-2) 

or  six  convolutions  make  up  the  last  term.   Since  each  of 
these  convolutions  is  of  the  type  previously  evaluated,  the 
solution  is  now  straightforward  and  the  last  term  becomes 

e^6  L  -1 

~A~  L     \  «^*->  ^*-\  ^v 

lO-l  u.^\  uvv  ' 
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_  XCjc>^cC^-L\ 


^-\ 


- — Wz# 


Z^-\ 


i^-l 


'prtf-   (.i^-Oy 


-  3       \   _J Cjos(zuvtyt   -  J (Lost    *  f  J 


to -i 


^o-^        3Ui- 


Lj\cofcfi^"J(  (c-5) 
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